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Abstract. - We study the direct enstrophy cascade in a two-dimensional flow generated 
in an electromagnetically driven thin layer of fluid. Due to the presence of bottom friction, 
the energy spectrum deviates from the classical Kraichnan prediction fc~ 3 . We find that the 
correction to the spectral slope depends on the thickness on the layer, in agreement with a 
theoretical prediction based on the analogy with passive scalar statistics. 



Introduction. Laboratory studies of two-dimensional turbulence are affected by the 
interaction with the three-dimensional environment in which they are embedded. This in- 
teraction is often frictional and generates an additional linear damping term to the equation 
of motion. A well known example, which is the object of the present letter, is the flow in a 
shallow layer of fluid, which is subject to the friction with the bottom of the tank where it is 
contained [1]. Another example is the motion of a soap film dumped from the interaction with 
the surrounding air [2]. Linear friction dumping is not restricted to laboratory experiments: 
an important example is the Ekman friction in geophysical fluid dynamics [3]. 

As predicted in a remarkable paper by Kraichnan in 1967 [4], when a two-dimensional 
layer of fluids is forced at intermediate scales, smaller than the domain-size and larger than 
the viscous dissipative length-scales, two different inertial ranges are observed. This is a 
consequence of the coupled conservation of energy (mean square velocity) and enstrophy (mean 
square vorticity) in the inviscid limit, a basic difference with respect to three-dimensional 
turbulence phenomenology. The energy flows toward the large scales giving rise to an inverse 
energy cascade characterized by the fc -5 / 3 Kolmogorov spectrum. Small scale statistics is 
governed by a direct enstrophy cascade which is expected to develop a smooth flow with fc -3 
energy spectrum, with a possible logarithmic correction [5] . 

The presence of a linear friction term affects the direct cascade in a dramatic way: the 
enstrophy flux across scales is no longer constant and the scaling exponent of the spectrum 
differs from the Kraichnan prediction fc~ 3 by a correction proportional to the friction intensity 
[6,7]. The origin of this correction can be understood exploiting the analogy between two- 
dimensional vorticity field and a scalar field passively advected by a smooth velocity field [8] . 
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Fig. 1 - a) A schematic upper view and side view of the experimental setup, b) Snapshot of the 
velocity field reconstructed by the FT technique on the regular at resolution 128 x 128 grid for the 
run 3 (see Table 0. 



For the case of a passive scalar with a finite lifetime it is possible to obtain explicit expressions 
for the scaling exponents of the power spectrum and structure functions which depend on the 
intensity of the linear dumping and on the statistics of finite-time Lyapunov exponents [9, 10]. 

In this letter we experimentally study the statistic of the direct enstrophy cascade in a 
two-dimensional flow generated in an electromagnetically driven thin layer of fluid. This 
experimental setup has been successfully used for studying several problems connected to 
two-dimensional turbulence [1, 11, 12]. Here we show that the friction with the bottom of the 
tank induces a correction on the spectral slope at large wave-numbers, in good agreement 
with theoretical predictions [6,7]. 

Experiment. - The experimental apparatus consists of a squared Plexiglas tank whose 
dimensions are 50 x 50 cm. 

The tank is filled with two layers of fluid, subsequentially injected from the bottom. The 
upper layer of fresh water is placed over a bottom layer of an electrolyte solution of water and 
NaCl (density ~ 1060<?/^). At the end of the filling procedure we obtain a stable stratification 
with controlled fluids thickness. In all the experiments the fluid thickness h\ for the lower fluid 
has been maintained equal to 0.3cm while the thickness of the upper fluid h,2 varies between 
0.3 and 0.7cm. The total thickness of the two layers will be denoted by h. 

Neodymium permanent magnets are placed below the bottom of the tank and disposed 
in four triangles with magnetic field in the vertical direction and alternating sign (see Fig^l. 
An electric current, horizontally driven through the cell, interacts with the magnetic field 
and moves the NaCl solution via the Lorenz force. The forcing current is generated by 
a computer controlled power supply which provides voltage signal of fixed amplitude and 
randomly alternating direction. The correlation time for voltage sign reversal is tipically 
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Fig. 2 - Decay of kinetic energy for the unforced runs (run 1 (+), run 2 (x), run 3 (o)). The lines 
represent exponential fits which give the friction coefficients reported in Table Q 



4s. The combined action of electric and magnetic forcing on the NaCl solution induces the 
continuous formation of opposite signed vortical structures whose characteristic length-scale 
Lf is related to the distance between opposite-signed magnets and whose characteristic time- 
scale is of the order of voltage sign reversal. 

The free fluid surface is seeded using tiny buoyant styrene particles (with typical size d ~ 
250/im), and the test section is illuminated using an array of lamps placed orthogonally to one 
side of the tank. The fluid flow is recorded using a standard speed video camera. A maximum 
duration of 6 minutes is chosen for each experiment, to ensure that both density stratification 
and two-dimensionality are maintained [13]. The velocity field has been reconstructed by 
image analysis based on a Feature Tracking (FT) approach [14]. This tracking procedure 
allows for higher seeding densities than classical Particle Tracking Velocimetry, and provides 
an accurate reconstruction of a large number (almost 20000) of Lagrangian trajectories. The 
interpolated Eulerian velocity field is therefore highly detailed, maximizing the information 
content of raw data. In Fig. ^ we show an example of instantaneous velocity field reconstructed 
at resolution 128 x 128. 

Theory and experimental results. The dynamics of a thin layer of fluid electromag- 
netically forced is described by the linear dumped 2D Navier-Stokes equations, which can be 



Table 1 - Parameters of the experiments. Friction coefficient a, root-mean-square velocity u r 
root-mean-square vorticity ui rms , spectral index correction £ 



Run # 


h(cm) 


ais- 1 ) 


Urms(cmS *) 




£ 


1 


0.8 


0.0371 


1.324 


0.747 


0.49 


2 


0.9 


0.0591 


1.330 


0.640 


0.78 


3 


1.0 


0.0686 


0.790 


0.602 


1.02 
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Fig. 3 - Kinetic energy spectrum for the run 1 (+), The dashed line represents the Kraichnan 
prediction fc -3 . In the inset we show the spectrum compensated with k~ 3 ~^ (solid line) with exponent 
£ given in Tableland with the Kraichnan prediction k~ 3 (dashed line). 



written for the vorticity lo — V x v as 

^ + v ■ Vuj = vV 2 uj - auj + f(t) (1) 

where v is the fluid viscosity, f(t) is the time-dependent external forcing. The bottom drag 
is parameterized by the linear friction term —au. Assuming a Poiseuille-like vertical velocity 
profile, the intensity of the friction coefficient can be related to the total thickness of the layer 
h as a oc v/h 2 [15]. 

In the inviscid-unforced limit (a = v = f = 0), the NS equations conserves both the 
kinetic energy E — l/2(v 2 ) and the enstrophy Z = 1/2(cj 2 ) [16]. For a > the enstrophy 
cascading to small scales is removed by friction, allowing to disregard the viscous term in 
(UJ [6]. In this limit, both the energy and the enstrophy decay exponentially in the unforced 
case with a decaying characteristic time r = 1/ (2a) 

E(t) = E(0)e~ 2at Z(t) = Z(0)e~ 2at (2) 

In Figure [21 we plot the decay of the total energy for three experiments with different total 
thickness h starting from time at which the electric forcing is switched off. The agreement with 
the exponential decay is remarkable and allow a direct measurement of the friction coefficient 
a. Further this is a confirmation a posteriori of the irrelevance of the viscous term. 

The remarkable prediction made in [7] is that for any a > the Kraichnan scaling exponent 
in the direct cascade has a correction proportional to a. The argument is based on an analogy 
with the dynamics of a scalar with finite lifetime t = a -1 passively transported by a smooth 
flow, which is governed by an equation formally identically to Q [19]. We remark that the 
extension of the passive scalar argument [10] to the active vorticity field is not trivial, as oj 
determines the velocity field. For completenes, in the following we report the "mean field" 
derivation of the exponent correction, a complete derivation for the active case can be found 
in [8]. 

Let us consider a fluctuation of vorticity Sui(Lf 7 0) = Q generated by the forcing at scale 
Lf and time 0, i.e. a blob of vorticity of size Lf at the rms value fi. Due to the chaotic 
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Fig. 4 - Kinetic energy spectrum for the run 2 (x). The dashed line represents the Kraichnan 
prediction fc -3 . In the inset we show the spectrum compensated with fc -3- ^ (solid line) with exponent 
£ given in Tableland with the Kraichnan prediction fc -3 (dashed line). 



flow, the blob is exponentially stretched with a mean rate given by the Lyapunov exponent 
A. Because of the incompressibility of the velocity field, after a time t the blob is contracted 
in the transverse dimension to a scale r = Lf exp(— Xt). This is the mechanism for the direct 
cascade phenomenon, i.e the fluctuation has been transported from the large injection scale 
Lf to the small scale r. Taking into account the decay induced by the friction we can write 



a/A 

Suj(r,t) ~ 5uj(L f ,Q)e- at ~Q ^-j (3) 

In the statistically steady regime sustained by the forcing, this argument provides the scaling 
exponent for the second-order vorticity structure function: 

5 2 (r)^((fcW) 2 )-« 2 (^) 2Q/A (4) 

and thus the scaling exponent of the enstrophy spectrum Z(k) ~ k~ x ~ 2a l x . Finally one obtain 
the mean field prediction for the energy spectrum [7] 

E{k) ~ fc~ 3 ^ (5) 

with correction exponent £ = 2a/ X. We remark that for a > the correction gives an energy 
spectrum steeper that k~ 3 which is a posteriori consistent with the assumption of smooth 
velocity field. 

A more refined version of this "mean-field" argument, which takes into account the fluctua- 
tions of the Lyapunov exponents [18], can be made [7,8]. The result is again a correction which 
makes the spectrum steeper than the Kraichnan prediction. An extensive numerical study 
of two-dimensional Navier-Stokes equations with friction has confirmed the above theoretical 
picture [8]. 

Figures |21 El and [5] shows the energy spectra obtained from Fourier transform of the FT 
velocity data for three experiments with different total thickness in stationary conditions. 
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Fig. 5 - Kinetic energy spectrum for the run 3 (o). The dashed line represents the Kraichnan 
prediction fc -3 . In the inset we show the spectrum compensated with fc -3- ^ (solid line) with exponent 
£ given in Tableland with the Kraichnan prediction fc -3 (dashed line). 

The spectra are obtained by averaging over about 50 realizations of the velocity field. A clear 
cascade range with power-law scaling is evident in intermediate wave-numbers 0.7cm -1 < k < 
4.0cm _1 . The forcing wavenumber kf ~ 0.7cm -1 corresponds to an injection scale Li ~ 9 cm, 
consistent with the size of the array of magnets. 

A fit of the energy spectra in the range 1.0cm -1 < k < 4.0cm -1 gives the exponent 
corrections £ ~ 0.49, £ ~ 0.78 and £ ~ 1.02 for the three runs respectively. A direct comparison 
with the theoretical prediction © would require the knowledge of the Lyapunov exponent of 
the flow. It is possible to give a simple estimation by considering the characteristic time in the 
direct cascade which is given by w~^ s . Therefore, we can assume that the Lyapunov exponent 
of the flow is proportional to uj rmsi a quantity which is easily determined from the velocity 
field. Finally, by considering a couple of different runs we have: 



Using the values of ctfj and Wj from table|]]we obtain the predictions £1/^3 = 0.81 and £2/^3 = 
0.44 which are close to the direct estimation from the spectra £1/^3 = 0.78 and £1/^3 = 0.49. 

Conclusion. - In this letter we have presented an experimental study of the effects of 
bottom friction on the direct enstrophy cascade observed in thin layer of fluid clectromagneti- 
cally forced. Direct measurements of the friction coefficient are obtained from the exponential 
decay of the total energy when the forcing is switched off. In the stationary forced case the 
energy spectra of the reconstructed velocity fields display a power law behavior, with a slope 
fc -3-, » which differs from the Kraichnan prediction fc -3 . The correction £ to the spectral slope 
is due to the friction exerted by the bottom wall on the fluid and increases with the friction 
intensity. Its value can be predicted by theoretical arguments, in good agreement with our 
experimental measurements. Being the correction roughly proportional to the square of the 
inverse of the total thickness of the fluid, the observed effect is expected to be extremely 
relevant in the case of thin layers. 
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